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In mean-field magnetohydrodynamics the mean electromotive force due to velocity and magnetic field fluctuations plays a 
crucial role. In general it consists of two parts, one independent of and another one proportional to the mean magnetic field. 
The first part may be nonzero only in the presence of mhd turbulence, maintained, e.g., by small-scale dynamo action. It 
corresponds to a battery, which lets a mean magnetic field grow from zero to a finite value. The second part, which covers, 
e.g., the a effect, is important for large-scale dynamos. Only a few examples of the aforementioned first part of mean 
electromotive force have been discussed so far. It is shown that a mean electromotive force proportional to the mean fluid 
velocity, but independent of the mean magnetic field, may occur in an originally homogeneous isotropic mhd turbulence if 
there are nonzero correlations of velocity and electric current fluctuations or, what is equivalent, of vorticity and magnetic 
field fluctuations. This goes beyond the Yoshizawa effect, which consists in the occurrence of mean electromotive forces 
proportional to the mean vorticity or to the angular velocity defining the Coriolis force in a rotating frame and depends on 
the cross-helicity defined by the velocity and magnetic field fluctuations. Contributions to the mean electromotive force 
due to inhomogeneity of the turbulence are also considered. Possible consequences of the above and related findings for 
the generation of magnetic fields in cosmic bodies are discussed. 
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Mean-field magnetohydrodynamics has proved to be a use- 
ful tool for studying the behavior of mean magnetic fields in 
turbulently moving electrically conducting fluids (see, e.g., 
Moffatt 1979, Krause & Radler 1980, Brandenburg & Sub- 
ramanian 2005). Within this framework both the magnetic 
field B and the fluid velocity U are split into mean parts, B 
and U, and fluctuating parts, b and u. Starting from the in- 
duction equation governing B it is concluded that the mean 
magnetic field B has to obey 

d t B = v V 2 B + V x(UxB + S), V -2* = 0.(1) 

Here, i] means the magnetic diffusivity of the fluid, for sim- 
plicity considered as independent of position, and £ the 
mean electromotive force caused by the velocity and mag- 
netic fluctuations, 



£ = {u x b) 



(2) 



Mean fields are defined by some kind of averaging satisfy- 
ing the Reynolds rules. They are denoted either by overbars 
or synonymously by angle brackets. 

The induction equation governing B also implies 



d t b = 77 V 2 6 



V x [(it x b)' 4 
V 



U x b 

6 = 0, 



u x B] 



(3) 



* Corresponding author: e-mail: khraedler@arcor.de 



where (uxb)' = uxb — (uxb). With this in mind we may 
conclude that £ can be represented as a sum of two parts, 



£ = £ W +£ (B) 



(4) 



where £^ is a functional of u and U, and £^ B ^ a functional 
of u, U and B, which is linear in B but vanishes if B is 
zero everywhere and at all past times (see, e.g., Radler 1976, 
2000, Radler & Rheinhardt 2007). These statements apply 
independently on whether or not u or U depend on B. If 
they depend on B and the total variation of £ with B is 
considered, £^ may well depend on B, and £^ need not 
be linear in B. 

A non-zero £^ corresponds to a battery. Assume for a 
moment that equation (1) for B with £ = has no grow- 
ing solutions. If then £' ' takes non-zero values, but £^ B ^ 
remains equal to zero, B grows, even if initially equal to 
zero, to a finite magnitude determined by £^°\ If, on the 
other hand, £^ remains equal to zero a non-zero £^ B ) may 
allow (if it has a suitable structure) a dynamo, that is, let 
an arbitrarily small seed magnetic field B grow exponen- 
tially (in the absence of back-reaction on the fluid motion 
even endlessly). A small non-zero £^ may deliver a seed 
field for such a dynamo. This possibility has been already 
discussed in the context of young galaxies (Brandenburg & 
Urpin 1998). 

In most of the general representations and applications 
of mean-field magnetohydrodynamics the part £^ ' of the 
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electromotive force £ has been ignored. Indeed, if it occurs 
at all, it decays to zero in the course of time except in cases 
in which an independent magnetohydrodynamic turbulence 
exists, e.g., as a result of a small-scale dynamo. 

The possibility of a non-zero £^ due to local, that is 
small-scale, dynamos in the solar convection zone has been 
discussed by Radler (1976). We express his statements here 
by 

£ (0) = c 7 7 + c n fi + c 7 n7 x (5) 

where 7 is a gradient, e.g., of the turbulence intensity, H the 
angular velocity responsible for the Coriolis force, and c 7 , 
cq and c 7 n some coefficients. More precisely, c 7 and c 7 n 
are scalars and cq is a pseudoscalar. 

Another interesting result has been derived by Yoshizawa 
(1990, see also Yoshizawa 1993 or Yoshizawa, Itoh & Itoh 
2003). Considering an originally homogeneous isotropic mag 
netohydrodynamic turbulence under the influence of a mean 
flow or a rigid-body rotation, or both, he found 

£ (0) = c w W + c f2 , (6) 



V B = 



(7) 



V -U = 0. 



where W = V x U is the mean vorticity, fi again the angu- 
lar velocity responsible for the Coriolis force, and cw and 
cq are pseudoscalars coefficients which are, roughly speak- 
ing, proportional to the cross-helicity (u ■ b). This result 
has recently been used for an interpretation of the Archontis 
dynamo (Sur & Brandenburg 2009). 

The main purpose of this paper is to demonstrate that 
a mean electromotive force £^ proportional to the mean 
fluid velocity U may occur in originally homogeneous isotropic 
magnetohydrodynamic turbulence. This should be expected 
as soon as there is a non-zero correlation between the fluc- 
tuating parts of velocity and electric current, u and j = 
yUg^V x b, or, what is equivalent, between the fluctuating 
parts of vorticity and magnetic field, uj = V x u and b; 
as usual, /j,o means the magnetic permeability. We express 
this condition roughly by saying that (u ■ j) or (w • 6) have 
to be unequal to zero. Unlike (u ■ b), which characterizes 
the linkage between vortex tubes and magnetic flux tubes, 
(u ■ j) quantifies the linkage between vortex tubes and cur- 
rent tubes. 

In Section 2 we explain the basis of our calculations and 
provide general relations for the determination of the mean 
electromotive force £^°\ In Section 3 we derive results for 
homogeneous isotropic turbulence, in particular the last- 
mentioned one, and we also reproduce that given by (6). 
Proceeding then in Section 4 to inhomogeneous turbulence 
we report on results related to those indicated in (5). The 
relevance of the results obtained in this paper and the need 
of further work are discussed in Section 5. 



frame. It is assumed that B and the fluid velocity U are 
governed by 

dtB = r]V 2 B + V x{UxB + H), 
d t U + U-VU = -g- 1 VP + vV 2 V 
-2H xU + F , 
where rj is again the magnetic diffusivity, v the kinematic 
viscosity and g the mass density of the fluid. P is the sum 
of hydrostatic and centrifugal pressure and f2 the angular 
velocity defining the Coriolis force. The external electromo- 
tive force H and the external ponderomotive force F will 
allow us to mimic magnetohydrodynamic turbulence. For 
the sake of simplicity we have ignored the back-reaction of 
the magnetic field on the fluid motion. 

Adopting the mean-field concept and taking averages of 
(7) we arrive at equations for the mean fields B and U. The 
equation for B differs from (1) only in so far as instead of 
U x B the sum U x B + H occurs. For the magnetic and 
velocity fluctuations b and u we further may derive 

d t b = r]V 2 b + V x [(u x b)' + U x b + u x B + h] 

V-b = 

d t u = -g~ 1 Vp + uV 2 u- (u- Vu)' (8) 
20 xu-U - Vu-u- VU + f 
V u = 0, 

where again (u x b)' stands for u x b (u x b), analogously 
(it • Vu)' for u ■ Vit — (u ■ Vu), and h, p and / are the 
fluctuating parts of H, P and F. The equations for b differ 
. from (3) only by the additional electromotive force h. 

We strive to calculate the part £^ of the mean electro- 
motive force. So we put simply B = in (8). This is not 
germane to the following considerations and could always 
be justified by choosing a suitable H. Basically the remain- 
der of equations (8) with a given U implies the possibility of 
a small-scale dynamo, that is of non-decaying b, even if h 
is equal to zero. In what follows we introduce however some 
further simplifying assumptions which undermine this pos- 
sibility, and we mimic a small-scale dynamo with a proper 
non-zero h. 

Let us assume that u and b depend only weakly on U 
and r2 so that £^ is linear in these quantities. We further 
assume that U varies only weakly in space and time so that 
£^ in a given point depend only on U and its first spatial 
derivatives in this point. Thus we have 



£, 



(o) 



£l m + a. 



^ipq^U p I dXq -\- CipQp . 



(9) 



,(00) 



as well as the coefficients a, p , b ipq and c ip 



where £\ 

independent of U and tt. 

We now split b and u according to 

b = + 6 (1) 



are 



+ ••• , U = U (0) +••• (10) 

into parts b^ and independent of U and H, parts b^ 
and of first order in U or fl and higher-order contri- 
butions, which are however not considered in what follows. 
The assumption of the linearity of £^ in U and l~i implies 



2 General concept 
2.1 Basic equations 

We consider a magnetic field B in a homogeneous incom- 
pressible electrically conducting turbulent fluid in a rotating 



= («(») x 



bW) + {«Wx6(°'). (11) 
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Returning to equations (8) and considering h and / as 
independent of U and fi, we find for b (0) and 

dtb<V = ?yW > + V x [(«(») x b^)' + h] 

V • &<°> = o 

d t U^ = -g-^ V P W + !/VV°> (12) 

-(/'■V/')'!/ 
V-uW =0. 
In the following we denote the turbulence defined by 

b (0) 

and it' ) as "background turbulence". In the equations re- 
sulting for b^ and we introduce some generalized 
second-order correlation approximation, that is, neglect all 
terms originating from (it x b)' and (it • Viz)'. Hence we 
have 

d t b (1) = v V 2 b (1) + V x (U x b (0) ) 

V ■ bW = 

<W X) = --Vp (1) + z/V 2 w (1) (13) 
Q 

-2(1 x u (°) - U ■ Vu(°) - tx (0) • Vt7 

V -u (1) =0. 

2.2 Relation for £ (0) 

In the following derivations we use a Fourier transformation 
of the form 

F(x,t) = JJ F(k,uj) cxp[i(k-x-ujt)]d 3 kduj. (14) 

The integrations are over all k and u>. 

In view of the determination of £^ we first note 

(u(x,t) x b(x,t))i = e ijk JJ Q jk (x,t-.k,uj)d 3 kduj , 

(15) 

where Qj k (x, t; k, uj) is the Fourier transform of 

Q jk (x,t;£,T) = 

(Uj (x + £/2, t + r/2) b k (x - |/2, t - r/2)) (16) 

with respect to £ and r. Adopting the formalism of Roberts 
& Soward (1975) we find that 

Q jk (x,t; k,oj) 

{ujik + k'/2, uj + J j2) b k (-k + k'/2, -u + uj' 12 ) 

exp[i(fc' x -oj't)]d 3 k'duj' ; (17) 

see Appendix A. As a consequence ofV-it = V- fo = 
the conditions 

VjQj k +2ikjQ jk = , V k Qj k - 2\k k Q jk = (18) 

have to be satisfied. Note that for the determination of £^ 
only the antisymmetric part of Qj k is needed. 

Considering £^ we restrict our attention simply to x = 
0. In that sense we put 



//« 



We consider the equations (12) as solved, that is, b^ and 
«( ) as known. Subjecting then the equations (13) for b^ 
and with U specified according to (19) to a Fourier 
transformation and eliminating the pressure term in the usual 
way we find 

i>V = -E[xk m U m bf ] - U im bW - k m U mn dbf } /dk n ] 
E = (T]k 2 - iujy 1 , biki=0 
= -N[ik m U m uf ] +U im uM 

-k m U mn (2hu ( ^/k 2 + duf ] /dk n ) (20) 

+2 (k-n)u^/k 2 ] 

N = (vk 2 - iujy 1 , Uiki = i). 
Calculating now Qj k on the basis of (17) and (20) we 
neglect again all contributions of higher than first order in 
U and ft. We further discard terms with more than one spa- 
tial derivative, in particular products of Uij with any other 
spatial derivative. Since Qj k should only weakly vary with 
x we expand (iijbk) in (17) for small k! and arrive so at 



Q jk = Q? k ) +i(E*-N) (k-U)Q^ 



+E*U km Qfl ~ NU jm QZ + 2NU mn k, J k m Q^> 
+\(E*' + N')U mn k m k n Q { $ Ik 2 
+ \{E* + N)U mn k m dQf k ] /dk n 
-2Ne imn k m (k-n)Q { ° k } /k 2 
-±(E*+N) (U-V)Q^ 



(o) 



l(0) 



(21) 



-\{E*' + n') (k-u) (fc.v)gSVfc 2 

+ie jmn [N((k ■ fi)V m + k m (n ■ V))Qf k /k 2 
-{2N - N')k m {k • 0)(fe ■ v)gi fc } /fc 4 ] . 

Here Qf k means Qj k for the background turbulence, that is, 

with u and b replaced by u/ ' and b^ ' . For simplicity argu- 
ments are dropped; those of Qj k and are (x, t; k, uj), 
those of E, N, etc. are (k, uj). The asterisk means complex 
conjugation and F' = kdF/dk. 

Returning now to the representation (9) of £^ we find 



£•(00) _ c 
^4 — fc i 



e«* J J Q% ] d 3 kdu; 
a lp = e ijk JJ [i(E* - N)k p - \{E* + N)\7 p 

-\{E*' + N'){k p /k 2 ){k- V)]Qf k ] d 3 kduj 
b ipq - e ijk JJ [E*5 kp Q { ^ + \{E* + N)k p 8Q jk /dk q 



-N(S JP - kjkp/k 2 )^] d 3 kduj 
■■ JJ [(N/k 2 ) (2hk p - i(h\7 p + k p V l )) Q 



(22) 



(0) 



(0) 



(19) 



+i(N/k 2 ){k p V k + k k V p )Q\ u k 
-i(2N- N')(k p /k 4 )(k ■ V) 

(k k Qf k ] -k t Qf)]d 3 kduj, 
where again the above remarks on arguments apply. 
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3 Homogeneous isotropic turbulence 
3.1 General result 

Consider now the simple case in which the background tur- 
bulence is homogeneous and isotropic and return first to 
(9). Since there is no isotropic vector we have £^ ^ = 0. 
Isotropy further implies a ip = cu5 ip , b^ k = c w eij k and 
Ci P = cnSij. Hence we obtain 



= Cu U + c w V x U + cq fi 



(23) 



with a scalar cjj and pseudoscalars cw and cq. 

Due to homogeneity and isotropy of the turbulence we 
have 



(5 jk 



k > kk -)^(k,u)-ie 3kl ^¥^(k,u) 



(24) 



where l>( ) and <&(°) are the Fourier transforms of 

$(°) = 

(ix (0) (x + £/2, t + r/2) ■ b (0) (x - £/2, t - r/2)> 
* (0) = (25) 
W(« (0) (x + €/2, « + r/2) • j<°> (as - €/2, t - r/2)) 

with respect to £ and r, and fioj^ stands for V x . The 
homogeneity implies that <j>(°) and \I/( ) as well as <I>( ) and 
are independent of x and further that 

Mo (« (0) (* + £/2, t + r/2) • j<°> (x - €/2, t - r/2)) = 
(w (0) (a: + €/2, t + r/2) • (x - |/2, t - r/2)) , (26) 

where w(°) is the vorticity of the velocity field «(°), that 
is w(°) = V x «(°). In general Q jk , $(°), *(°) as well 
as $(°^ and ^(°) may depend on t. If we however assume 
that the turbulence shows in addition to its homogeneity also 
statistical steadiness this dependence vanishes. In addition 
the arguments (x + £/2, t + r/2) and (x - £/2, t - r/2) 
in (25) and (26) may then be replaced, e.g., by (x,t) and 
(x — £, t — r) or by (x + £, t + r) and (x, t), respectively. 
With (22) and (24) we find 

cu = (E* - N) *W kd 3 kdcu 

c w = ^JJ (2E* + N) $(°) d 3 kdu 
N& 0) d 3 kduj. 



(27) 



We point out that E = (27r) 4 GW and iV = (2tt) 4 GH 
where the G (7) are the well-known Green's functions de- 
fined by 

GW^.t) = (47r 7 r)- 3 / 2 exp(-e 2 /47r) for r > 
G (7) (C,^) = forr < 0. (28) 



Considering this and applying the convolution theorem to 
(27) we obtain 



(u 



(0), 



(x,t)-j« > \x-t,t-T))d !l ZdT 



(« (0) (x, t) • b (0) (a; - £, t - r)) d 3 £dr 

(u<°> (a;, t) • 6 (0) (x - €, t - r)) d 3 ^ dr . 
Here the integrations are over all £ and primarily also over 
all r. However, since the G {r)) = G {u) = for r < 0, they 
involve in fact only positive r. 

The most remarkable result of our derivations is that a 
contribution to £^ proportional to U, that is, a term cjjU 
in (23), may occur. This possibility has not previously been 
considered in the literature. According to (29), as long as 
(u/°) (x, t)-j^ (x— £, t—r)) does not vanish and 77 ^ j/,the 
coefficient cjj may well be different from zero. In the special 
case 77 = v, however, it is equal to zero. In what follows the 
occurrence of that contribution to £^ proportional to U is 
labeled as "(u ■ j) effect". 

For non-vanishing (u^ (x,t) ■ (x — £, t — r)) the 
pseudoscalar coefficients cw and cq will in general be dif- 
ferent from zero. Then, as already found by Yoshizawa (1990), 
contributions to £^ proportional to the mean vorticity Vx[/ 
and to the angular velocity f2 will occur. We refer to them 
as "(it • b) effects" or "Yoshizawa effects". 

We stress that the (u ■ j) effect is well possible under 
circumstances in which (u^(x, t) ■ b^°\x — £, t — r)) is 
equal to zero, that is, in which there are no (u ■ b) effects. 

3.2 Special cases 

Let us first consider £^ in some limiting cases with respect 
to r/ and v. We use the fact that 

G ( %,r) ^ 3 (£) as 7^0. (30) 

In the limit defined by 77 — > and c^oowe obtain 

1 , !^ 
Cu = --A, c w — -C , c o = 0, 

in the limit 77 — > 00 and v — > 

112 

C(7 = X^' c w = a C '' co^-oC' ( 31 ) 

and in the limit 77, — > 



1 2 
c w = 0, c w = -G, cn = --C, (32) 



where 



A = Mo / (M (o) (^,0-J (O) (^^-T))dr 
Jo 

/•OO 

G = / (tt (0) (x,t) -b (0) (a;,t-r)))dr. 
Jo 



(33) 
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Instead of the last relations we may also write 

A = » (uW(x,t)-jW(x,t))T A 
C = (u^(x,t)-b^(x,t)))T C 



(34) 



with correlation times ta and tq defined by equating the 
respective right-hand sides of (33) and (34). 

In view of a numerical test we also consider the case 
in which («(°) (x,t) ■ j (0) (x-£,t- r)> and («(°) (aj, t) ■ 

(x — £, t — t)), as far as they enter into the integrals in 
(29), do not markedly vary with r. Ignoring the dependence 
on r completely and using 



f°° 1 



(35) 



we find 



c -^G + ^) ct ' c *=-<L c " (36) 



with 



A = [ ( M (°)(i,t).(Vxb<°)(x-^)))^ 

J oo S 

Ct = / (tt(°)(x,t).6(°)(x-€,t))^. (37) 

J oo S 

We may introduce vector potentials and such that 

V x ^(°) = «(°) , V • ^<°> = 

V x a <°> = b (0) , V • o(°) = 

and therefore 

a (o) {x) = J_ I Vxb (o) (a; _^^i. 



(38) 



(39) 



For simplicity the argument t is omitted everywhere. The 
result (36) can then be written in the form 



cu = 



c w 





1 




1 


1 




3 


( 


V 


V 


1 










3 


( 


V 


+ 


2v , 



A* 



cn = -— (40) 
2v 



with 



A* = (tt(°) • aW) = (V (0) • b (0) ) 

C-t = (^(°).6( )) = ((Vx U (°)).a( )). (41) 

The arguments of the quantities in the angle brackets are, of 
course, always (x,t). 



3.3 A numerical test 

As a check of the above derivations, the electromotive force 
£^ and so the coefficient cu have been determined with 
numerical solutions of the equations (8). For these calcu- 
lations B has been put equal to zero, and U was spec- 
ified via the initial condition to be constant in space and 
turned out to remain nearly constant in time, too. The forc- 
ing fields h and / were taken as periodic in the space co- 
ordinates x, y and z, and steady. More precisely, h and 
/ differed only by constant factors from the vector field 
e(kx) = (sin kz, sinfca;, sin ky), with a constant k. The 
flow which would result from / is the no-cosine ABC flow 
of Archontis (2000, see also Dorch & Archontis 2004 and 
Cameron & Galloway 2006). Flows of this type are non- 
helical. For this reason the (u ■ j) effect should occur, but 
no (u ■ b) effects. Corresponding to the steadiness of h and 
/ only steady b and u were considered. The average which 
defines mean fields was taken over all x, y and z or, equiva- 
lent to this, over a periodic box. No approximation such as, 
e.g., the second-order correlation approximation was used. 

Let us specify the result for cjj given by (40) and (41), 
which has been derived in the second-order correlation ap- 
proximation, to the described situation. Relying on (12) we 
assume that and b^ are dominated by contributions 
proportional to e(kx). So we find 



cu = c R m 1 - 



1 



co = -- 



«o (« (0) -j (0) ) 



(42) 

with the magnetic Reynolds number R m and the magnetic 
Prandtl number P m defined by 



(43) 



and ui°2 s 



u (0)2)l/2_ 

In Table 1 and Figure 1 the numerically determined val- 
ues of cu/cq are given in dependence of R m and P m . In 
agreement with what we have found in our analytical calcu- 
lations in the second-order correlation approximation their 
signs change with growing P m at P m = 1. Moreover, in 
most cases the numerically determined values completely 
agree with those obtained in this approximation, that is, with 
i? m (l — 1/ P m ). Deviations occur only if the fluid Reynolds 
number _R C = u^ s /vk, that is R c = R m /Pm, exceeds a 
value of about 5. We should emphasize that all our numer- 
ical solutions are laminar and perfectly regular, just as in 
Figure 2 (upper row) of Sur & Brandenburg (2009). 

4 Inhomogeneous turbulence 

Let us add some results for the case in which the turbulence 
is no longer homogeneous and so also no longer isotropic. 
We admit now that quantities like (w (0) (a; + £/2, t + r/2) • 
j (0) (;E-£/2,t-T/2))and (-it(°)(a;+£/2,t+T/2)-b (0) (;E- 
£/2, i — r/2)} may depend on x so that their gradients with 
respect to x do not generally vanish. Then, in addition to 
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Table 1 Numerically calculated values of cu /co for several P m 
and P m , to be compared with the values derived in the second- 
order correlation approximation, P m (l — 1/P m ). 



Rm 




cu/co 


P m (l - 1/P m ) 


0.2 


4 


0.15 


0.15 




2 


0.10 


0.10 




1 


2.6 ■ 10~ 6 







0.2 


-0.80 


-0.80 




0.04 


-4.80 


-4.80 




0.01 


-15.5 


-19.8 


1 


5 


0.80 


0.80 




1 


1.2 • 10 -6 







0.2 
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Fig. 1 Numerically obtained values of \cu\/co, indicated by o 
for P m < 1 (where cu < 0) but by • for P m > 1 (where cu > 0), 
and curves representing P m |l — 1/P m | for R m — 10 (solid), 
R m = 1 (dotted), and R m = 0.2 (dashed). 



the contributions to £^ given in (23), other contributions 
are possible. Symmetry considerations suggest 

f (o) 



cu U + cw V x U + cn fi 
+9 + 9u x U + g n x fi , 



(44) 



where cu, Cw and cq may now vary in space, g and g n are 
gradients of scalars of the type of cu, and g v is a gradient 
of a pseudoscalar of type of cw or cn. For the sake of sim- 
plicity only terms up to the first order in spatial derivatives 
are regarded here. 

We may determine the contributions to £^ mentioned 
in (44) again on the basis of (15) and (21). However, relation 
(24) for , which applies to homogeneous isotropic tur- 
bulence only, has to be modified. Admitting now that $(°) 



and * (0) defined by (25) and so also $ (0) and # (0) depend 
on x and assuming that their gradients are small we add 
terms which are linear in these gradients to the right-hand 
side of (24). Considering further the conditions (18) we find 



{5jk 
1 



e jkl ki{2\ + — (fe- V)) 
1 



(45) 



--^(kjCkim + k k e jlm )kiV m # (0) (a?; k,uj) . 

As in the case of homogeneous turbulence Qjk, and 
may depend on t, what is however of minor importance 
in this context and therefore not explicitly indicated. 

A straightforward calculation confirms then (44). As for 
Cu, cw and cn we find, as expected, again the relations (27), 
now with ^( ) and <|( ) depending in general on x, and so 
also (29). Furthermore the calculation yields 

9 = - x -v 

y 6 



*(°)fc- 2 d 3 fcdo; 
9U = l V jj(E* -N)¥ >d 3 kdu; 



(46) 



9n = -3V 



JV#(°) k- 2 d 3 kduj. 



Using again the convolution theorem in combination with 
above-mentioned connection between N and £?(") and be- 
tween E and G^) we arrive at the equivalent relations 



j (a \x-H/2,t-T/2))^dT 



9U = \V 



jj '(<?("> (£,t)-G(">(€,t)) (47) 
(w (0) {x + £/2, t + r/2) • 6 (0) (a; - £/2, t - r/2)} d 3 £dr 



Mo, 

67T 



(t*W (a; + |/2, t + r/2) • jW (x - £/2, t - r/2)) d d £ dr . 

These results confirm in some sense the statements made 
in the paper by Radler (1976), formulated above in (5). They 
show however that the vectors c 7 7 and c 7 n7 should not, 
as suggested there, be understood in the sense of V(ix 2 ). 
These vectors rather correspond to g or g n as given in (46) 
and (47). Roughly speaking, they should be interpreted in 
terms of V(u • j). 



5 Discussion 

The most remarkable result of our calculations is that the 
mean electromotive force £^ in a homogeneous isotropic 
magnetohydrodynamic turbulence may have a contribution 
proportional to the mean fluid velocity U, that is, = 
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cjj U + ■ ■ ■ . We have labeled the occurrence of this con- 
tribution as (u ■ j) effect. The coefficient cjj turned out to 
be in general unequal to zero if only a non-zero correlation 
exists between the fluctuating parts of the fluid velocity and 
the electric current density, u and j = /i^ 1 V x b, or be- 
tween the fluctuating parts of the vorticity and the magnetic 
field, u> — V x u and b. As far as the second-order correla- 
tion approximation applies, c v vanishes for r\ = u, that is, it 
changes its sign ifv/r] varies and passes through v/r] = 1. 

The occurrence of magnetohydrodynamic turbulence does 
not automatically imply non-zero correlations of u and j, 
or u> and b. It depends on the special circumstances whether, 
e.g., (u ■ j) or (u> ■ b) are different from zero and what their 
signs are. For this and other reasons more work is needed 
to explore the importance of the (u ■ j) effect in specific 
settings. In general the (u ■ j) effect is accompanied by the 
(u ■ b) effects. It has then also to be investigated which of 
these effects dominates. 

At first glance there seems to be a inconsistency of our 
result concerning the (u ■ j) effect in so far as £ (0) should 
not depend on the choice of the frame of reference but U ob- 
viously does. We must however keep in mind that we have 
fixed the frame of reference in our calculation by assuming 
that there isotropic turbulence occurs in the limit U — > 0. 
When estimating the (u -j) effect we have therefore to spec- 
ify U as the mean velocity of the fluid relative to the frame 
in which the assumed causes of turbulence (in simulations 
the forcing) would, in the absence of mean motion, just lead 
to isotropic turbulence. 

There is still another issue which has to be considered 
when applying our result concerning the (u ■ j) effect to 
a specific situation. The deviation of the turbulence from 
isotropy due to the homogeneous part of the mean motion 
is crucial for that effect. It has to be scrutinized whether 
such a deviation indeed occurs under the considered cir- 
cumstances. The sometimes assumed "Galilean invariance" 
of the turbulence (e.g., Sridhar & Subramanian 2009), that 
is, its independence of that part of the mean motion, would 
exclude the (u ■ j) effect. 

In principle we could have determined £^ in a frame 
in which U vanishes. An anisotropy of the turbulence, as 
calculated in the frame used above, would lead to a non- 
vanishing contribution instead of cjjU, which then 
has to be considered as another description of the (u ■ j) 
effect. 
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A Derivation of relation (17) for Q jk 

Start from Qjk{x, t; £, r) as given in (16) and introduce there the 
Fourier representations of Uj and b k so that 

Qjk(x,t;£,r) = 

// //<^'( fc W)MfcW)> 

exp (i((fe t + fe*) • x + (fe + - fe*) • |/2 (Al) 
-(u; f +w*)t- (w f - J>/2)) 

d 3 *^ dw + d 3 fe*du; t . 
Change then the integration variables according to 
tf=k + k'/2, k x = -k + k'/2 
a/ = oj + u//2 , uj % = -uj + uj' /2 (A2) 
and find so 

Q 3 k(x,t;$,r) = 

(uj(k + k'/2,ui + tj'/2) 

b k (-k + k' /2, -lu + cu' /2)) 

exp (i(fc' • x — u)'t)) d 3 fc' du/ (A3) 
cxp (i(fc ■ £ — ojt)) d 3 fc du; . 

This shows that Qjk given by (17) is indeed the Fourier transform 
of Qjk(x, t; £, r) with respect to | and r. 
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